Detecting episodic evolution through Bayesian inference of molecular clock models (preprint)

Molecular evolutionary rate variation is a key aspect of the evolution of many organisms that can be modelled using molecular clock models. For example, fixed local clocks revealed the role of episodic evolution in the emergence of SARS-CoV-2 variants of concern. Like all statistical models, however, the reliability of such inferences is contingent on an assessment of statistical evidence. We present a novel Bayesian phylogenetic approach for detecting episodic evolution. It consists of computing Bayes factors, as the ratio of posterior and prior odds of evolutionary rate increases, effectively quantifying support for the effect size. We conducted an extensive simulation study to illustrate the power of this method and benchmarked it to formal model comparison of a range of molecular clock models using (log) marginal likelihood estimation, and to inference under a random local clock model. Quantifying support for the effect size has higher sensitivity than formal model testing and is straight-forward to compute, because it only needs samples from the posterior and prior distribution. However formal model testing has the advantage of accommodating a wide range molecular clock models. In contrast, the random local clock had low power for detecting episodic evolution. In an empirical analysis of a data set of SARS-CoV-2 genomes, we find 'very strong' evidence for episodic evolution. Our results provide guidelines and practical methods for Bayesian detection of episodic evolution, as well as avenues for further research into this phenomenon.

Phytest: Quality Control for Phylogenetic Analyses (preprint)

Motivation: The ability to automatically conduct quality control checks on phylogenetic analyses is becoming more important with the increase of genetic sequencing and use of real-time pipelines e.g. in the SARS-CoV-2 era. Implementations of ”nowcasting” or real-time phylogenetic analyses require automated testing to make sure that problems in the data are caught automatically within analysis pipelines and in a timely manner. Here we present Phytest (version 1.0) a tool for automating quality control checks on sequence, trees and metadata during phylogenetic analyses. Results: Phytest is a phylogenetic analysis testing program that easily integrates into existing phylogenetic pipelines. We demonstrate the utility of Phytest with real-world examples. Availability: Phytest source code available on GitHub (https://github.com/phytest-devs/phytest) and can be installed via PyPI with the command ‘pip install phytest‘. Extensive documentation can be found at https://phytest-devs.github.io/phytest/. Contact: wytamma.wirth@unimelb.edu.au

Assessment of COVID-19 intervention strategies in the Nordic countries using genomic epidemiology (preprint)

The Nordic countries, defined here as Norway, Sweden, Denmark, Finland and Iceland, are known for their comparable demographics and political systems. Since these countries implemented different COVID-19 intervention strategies, they provide a natural laboratory for examining how COVID-19 policies and mitigation strategies affected the propagation, evolution and spread of the SARS-CoV-2 virus. We explored how the duration, the size and number of transmission clusters, defined as country-specific monophyletic groups in a SARS-CoV-2 phylogenetic tree, differed between the Nordic countries. We found that Sweden had the largest number of COVID-19 transmission clusters followed by Denmark, Norway, Finland and Iceland. Moreover, Sweden and Denmark had the largest, and most enduring, transmission clusters followed by Norway, Finland and Iceland. In addition, there was a significant positive association between transmission cluster size and duration, suggesting that the size of transmission clusters could be reduced by rapid and effective contact tracing. Thus, these data indicate that to reduce the general burden of COVID-19 there should be a focus on limiting dense gatherings and their subsequent contacts to keep the number, size and duration of transmission clusters to a minimum. Our results further suggest that although geographical connectivity, population density and openness influence the spread and the size of SARS-CoV-2 transmission clusters, country-specific intervention strategies had the largest single impact.

The emergence of SARS-CoV-2 variants of concern is driven by acceleration of the evolutionary rate (preprint)

The ongoing SARS-CoV-2 pandemic has seen an unprecedented amount of rapidly generated genome data. These data have revealed the emergence of lineages with mutations associated to transmissibility and antigenicity, known as variants of concern (VOCs). A striking aspect of VOCs is that many of them involve an unusually large number of defining mutations. Current phylogenetic estimates of the evolutionary rate of SARS-CoV-2 suggest that its genome accrues around 2 mutations per month. However, VOCs can have around 15 defining mutations and it is hypothesised that they emerged over the course of a few months, implying that they must have evolved faster for a period of time. We analysed genome sequence data from the GISAID database to assess whether the emergence of VOCs can be attributed to changes in the evolutionary rate of the virus and whether this pattern can be detected at a phylogenetic level using genome data. We fit a range of molecular clock models and assessed their statistical fit. Our analyses indicate that the emergence of VOCs is driven by an episodic increase in the evolutionary rate of around 4-fold the background phylogenetic rate estimate that may have lasted several weeks or months. These results underscore the importance of monitoring the molecular evolution of the virus as a means of understanding the circumstances under which VOCs may emerge.

Genomics-Informed Responses in the Elimination of COVID-19 in Australia (preprint)

BACKGROUND: A cornerstone of Australia’s ability to control COVID-19 has been effective border control, using an extensive supervised quarantine program. However, a rapid recrudescence in COVID-19 cases was observed in the state of Victoria in June 2020. Here, we describe the genomic findings that located the source of this second wave as a breach in supervised hotel quarantine and demonstrate the successful elimination of COVID-19 for a second time in Australia.METHODS: Genome sequencing was performed on all available SARS-CoV-2-positive samples in Victoria and integrated genomic and epidemiological investigation undertaken.RESULTS: At 31st January 2021, 20,451 COVID-19 cases were reported in Victoria; samples were sequenced from 75% of cases (15,431/20,451). Genomics revealed 98% (10,426/10,646) of locally-acquired cases during the second wave were derived from a single incursion from hotel quarantine, with the outbreak strain rapidly detected in other Australian states and territories. Phylodynamic analyses indicated an epidemic growth rate comparable to emerging variants, such as B.1.1.7 in the United Kingdom. Strict public health interventions resulted in the elimination of the outbreak strain by 29th October 2020. Subsequent cases represented independent international or interstate introductions, with limited local spread.CONCLUSIONS: Rapid escalation of clonal outbreaks can occur from even a single breach of control practices, as revealed through our genomic ‘enhanced outbreak-detection' system. The subsequent elimination and rapid control of new SARS-CoV-2 incursions reinforce that decisive public health responses to emergent cases are effective even with high epidemic growth rates, and “elimination” should be favored in settings where this is achievable.FUNDING STATEMENT: The Microbiological Diagnostic Unit Public Health Laboratory (MDU PHL) and the Victorian Infectious Diseases Reference Laboratory (VIDRL) at The Doherty Institute are funded by the Victorian Government. This work was supported by the National Health and Medical Research Council, Australia (NHMRC); Partnership Grant (APP1149991), Investigator Grant to BPH (APP1196103), Investigator Grant to DAW (APP1174555), Research Fellowship to TPS (APP1105525), MRFF COVID-19 Genomics Grant (MRF9200006).DECLARATION OF INTERESTS: None to declare. ETHICS APPROVAL STATEMENT: Data were collected in accordance with the Victorian Public Health and Wellbeing Act 2008. Ethical approval was received from the University of Melbourne Human Research Ethics Committee (study number 1954615.3).

The first wave of the Spanish COVID-19 epidemic was associated with early introductions and fast spread of a dominating genetic variant (preprint)

The COVID-19 pandemic has shaken the world since the beginning of 2020. Spain is among the European countries with the highest incidence of the disease during the first pandemic wave. We established a multidisciplinar consortium to monitor and study the evolution of the epidemic, with the aim of contributing to decision making and stopping rapid spreading across the country. We present the results for 2170 sequences from the first wave of the SARS-Cov-2 epidemic in Spain and representing 12% of diagnosed cases until 14th March. This effort allows us to document at least 500 initial introductions, between early February-March from multiple international sources. Importantly, we document the early raise of two dominant genetic variants in Spain (Spanish Epidemic Clades), named SEC7 and SEC8, likely amplified by superspreading events. In sharp contrast to other non-Asian countries those two variants were closely related to the initial variants of SARS-CoV-2 described in Asia and represented 40% of the genome sequences analyzed. The two dominant SECs were widely spread across the country compared to other genetic variants with SEC8 reaching a 60% prevalence just before the lockdown. Employing Bayesian phylodynamic analysis, we inferred a reduction in the effective reproductive number of these two SECs from around 2.5 to below 0.5 after the implementation of strict public-health interventions in mid March. The effects of lockdown on the genetic variants of the virus are reflected in the general replacement of preexisting SECs by a new variant at the beginning of the summer season. Our results reveal a significant difference in the genetic makeup of the epidemic in Spain and support the effectiveness of lockdown measures in controlling virus spread even for the most successful genetic variants. Finally, earlier control of SEC7 and particularly SEC8 might have reduced the incidence and impact of COVID-19 in our country.

The impact of early public health interventions on SARS-CoV-2 transmission and evolution (preprint)

Many countries have attempted to control COVID-19 through the implementation of non-pharmaceutical interventions. However, it remains unclear how different control strategies have impacted SARS-CoV-2 virus transmission dynamics at the local level. Using complete SARS-CoV-2 genomes, we inferred the relative frequencies of virus importation and exportation, as well as virus transmission chain dynamics in Nordic countries - Denmark, Finland, Iceland, Norway and Sweden - during the first months of the pandemic. Our analyses revealed that Sweden experienced more numerous transmission chains, which tended to have more cases, and were of longer duration, a set of features that increased with time. Together with Denmark, Sweden was also a net exporter of SARS-CoV-2. Hence, Sweden effectively constituted an epidemiological and evolutionary refugia that enabled the virus to maintain active transmission and spread to other geographic localities. This analysis highlights the utility of genomic surveillance where active transmission chain monitoring is a key metric.

Genomic epidemiology reveals transmission patterns and dynamics of SARS-CoV-2 in Aotearoa New Zealand (preprint)

New Zealand, a geographically remote Pacific island with easily sealable borders, implemented a nation-wide lockdown of all non-essential services to curb the spread of COVID-19. New Zealand has now effectively eliminated the virus, with low numbers of new cases limited to new arrivals in managed quarantine facilities at the border. Here, we generated 649 SARS-CoV-2 genome sequences from infected patients in New Zealand with samples collected between 26 February and 22 May 2020, representing 56% of all confirmed cases in this time period. Despite its remoteness, the viruses imported into New Zealand represented nearly all of the genomic diversity sequenced from the global virus population. The proportion of D614G variants in the virus spike protein increased over time due to an increase in their importation frequency, rather than selection within New Zealand. These data also helped to quantify the effectiveness of public health interventions. For example, the effective reproductive number, Re, of New Zealands largest cluster decreased from 7 to 0.2 within the first week of lockdown. Similarly, only 19% of virus introductions into New Zealand resulted in a transmission lineage of more than one additional case. Most of the cases that resulted in a transmission lineage originated from North America, rather than from Asia where the virus first emerged or from the nearest geographical neighbour, Australia. Genomic data also helped link more infections to a major transmission cluster than through epidemiological data alone, providing probable sources of infections for cases in which the source was unclear. Overall, these results demonstrate the utility of genomic pathogen surveillance to inform public health and disease mitigation.

Tracking the COVID-19 pandemic in Australia using genomics (preprint)

BACKGROUND: Whole-genome sequencing of pathogens can improve resolution of outbreak clusters and define possible transmission networks. We applied high-throughput genome sequencing of SARS-CoV-2 to 75% of cases in the State of Victoria (population 6.24 million) in Australia. METHODS: Cases of SARS-CoV-2 infection were detected through active case finding and contact tracing. A dedicated SARS-CoV-2 multidisciplinary genomic response team was formed to enable rapid integration of epidemiological and genomic data. Phylodynamic analysis was performed to assess the putative impact of social restrictions. RESULTS: Between 25 January and 14 April 2020, 1,333 COVID-19 cases were reported in Victoria, with a peak in late March. After applying internal quality control parameters, 903 samples were included in genomic analyses. Sequenced samples from Australia were representative of the global diversity of SARS-CoV-2, consistent with epidemiological findings of multiple importations and limited onward transmission. In total, 76 distinct genomic clusters were identified; these included large clusters associated with social venues, healthcare facilities and cruise ships. Sequencing of sequential samples from 98 patients revealed minimal intra-patient SARS-CoV-2 genomic diversity. Phylodynamic modelling indicated a significant reduction in the effective viral reproductive number (Re) from 1.63 to 0.48 after the implementation of travel restrictions and population-level physical distancing. CONCLUSIONS: Our data provide a comprehensive framework for the use of SARS-CoV-2 genomics in public health responses. The application of genomics to rapidly identify SARS-CoV-2 transmission chains will become critically important as social restrictions ease globally. Public health responses to emergent cases must be swift, highly focused and effective.

Temporal signal and the phylodynamic threshold of SARS-CoV-2 (preprint)

The ongoing SARS-CoV-2 outbreak marks the first time that large amounts of genome sequence data have been generated and made publicly available in near real-time. Early analyses of these data revealed low sequence variation, a finding that is consistent with a recently emerging outbreak, but which raises the question of whether such data are sufficiently informative for phylogenetic inferences of evolutionary rates and time scales. The phylodynamic threshold is a key concept that refers to the point in time at which sufficient molecular evolutionary change has accumulated in available genome samples to obtain robust phylodynamic estimates. For example, before the phylodynamic threshold is reached, genomic variation is so low that even large amounts of genome sequences may be insufficient to estimate the viruss evolutionary rate and the time scale of an outbreak. We collected genome sequences of SARS-CoV-2 from public databases at 8 different points in time and conducted a range of tests of temporal signal to determine if and when the phylodynamic threshold was reached, and the range of inferences that could be reliably drawn from these data. Our results indicate that by February 2nd 2020, estimates of evolutionary rates and time scales had become possible. Analyses of subsequent data sets, that included between 47 to 122 genomes, converged at an evolutionary rate of about 1.1x10-3 subs/site/year and a time of origin of around late November 2019. Our study provides guidelines to assess the phylodynamic threshold and demonstrates that establishing this threshold constitutes a fundamental step for understanding the power and limitations of early data in outbreak genome surveillance.

Direct RNA sequencing and early evolution of SARS-CoV-2 (preprint)

Fundamental aspects of SARS-CoV-2 biology remain to be described, having the potential to provide insight to the response effort for this high-priority pathogen. Here we describe the first native RNA sequence of SARS-CoV-2, detailing the coronaviral transcriptome and epitranscriptome, and share these data publicly. A data-driven inference of viral genetic features and evolutionary rate is also made. The rapid sharing of sequence information throughout the SARS-CoV-2 pandemic represents an inflection point for public health and genomic epidemiology, providing early insights into the biology and evolution of this emerging pathogen.